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Abstract. The gamma-ray astronomical observatory INTEGRAL, succesfully launched on 17 th October 2002, 
carries two large gamma-ray telescopes. One of them is the coded-mask imaging gamma-ray telescope onboard 
the INTEGRAL satellite (IBIS) which provides high-resolution (« 12') sky images of 29° x 29° in the energy 
range from 15 keV to 10 MeV with typical on-axis sensitivity of ~ 1 mCrab at 100 keV (3<j, 10 6 s exposure). 
We report here the general description of the IBIS coded-mask imaging system and of the standard IBIS science 
data analysis procedures. These procedures reconstruct, clean and combine IBIS sky images providing at the 
same time detection, identification and preliminary analysis of point-like sources present in the field. Spectral 
extraction has also been implemented and is based on simultaneous fitting of source and background shadowgram 
models to detector images. The procedures are illustrated using some of the IBIS data collected during the inflight 
calibrations and present performance is discussed. The analysis programs described here have been integrated as 
instrument specific software in the Integral Science Data Center (ISDC) analysis software packages currently used 
for the Quick Look, Standard and Off-line Scientific Analysis. 
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1. The IBIS coded aperture imaging system 

The IBIS telescope (Imager on Board of the INTEGRAL 
Satellite) (Ubertini et al. 2003), launched onboard the 
ESA gamma-ray space mission INTEGRAL (Winkler et 
al. 2003) on October 2002, is a hard-X ray/soft 7-ray 
telescope based on a coded aperture imaging system 
(Goldwurm et al. 2001). 

In coded aperture telescopes (Fig. 1) (Caroli et al. 
1987, Goldwurm 1995, Skinner 1995) the source radiation 
is spatially modulated by a mask of opaque and transpar- 
ent elements before being recorded by a position sensitive 
detector, allowing simultaneous measurement of source 
plus background (detector area corresponding to the mask 
holes) and background fluxes (detector corresponding to 
the opaque elements). Mask patterns are designed to al- 
low each source in the field of view (FOV) to cast a unique 
shadowgram on the detector, in order to avoid ambiguities 
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in the reconstruction of the sky image. This reconstruction 
(deconvolution) is generally based on a correlation proce- 
dure between the recorded image and a decoding array 
derived from the mask pattern. Assuming a perfect posi- 
tion sensitive detector plane (infinite spatial resolution), 
the angular resolution of such a system is then defined 
by the angle subtended by one hole at the detector. The 
sensitive area instead depends on the number of all trans- 
parent elements of the mask viewed by the detector. In the 
gamma-ray domain where the count rate is dominated by 
the background, the optimum transparent fraction is one 
half. The field of view (sky region where source radiation 
is modulated by the mask) is determined by the mask and 
the detector dimensions and their respective distance. To 
optimize the sensitive area of the detector and have large 
FOVs, masks larger than the detector plane are usually 
employed. The FOV is thus divided in two parts. The fully 
coded (FC) FOV for which all source radiation directed 
towards the detector plane is modulated by the mask and 
the Partially Coded (PC) FOV for which only a fraction 
of it is modulated by the mask. The rest, if detected, can- 
not be easily distinguished from the background. If holes 
are uniformly distributed the sensitivity is approximately 
constant in the FCFOV and decreases in the PCFOV. 

Representing the mask with an array M of 1 (transpar- 
ent) and (opaque) elements, the detector array D is given 



Goldwurm et al.: IBIS Science Data Analysis 



\ \ 



G 




Fig. 1. Scheme (not in scale) of the IBIS imaging sys- 
tem: a coded mask of 11.2 x 11.2 x 16 mm 3 tungsten 
elements, an upper low-energy detector plane (red) of 
128x128 CdTe crystals (4x4x2 mm 3 , pitch pixel-to- 
pixel = 4.6 mm), a high energy bottom detector plane 
(green) of 64x64 Csl bars (8.4 x 8.4 x 30 mm 3 , pitch 
pixel-to-pixel = 9.2 mm), an active BGO veto system 
(blue) composed of 19 crystals viewed by photomulti- 
pliers, a passive shield made by a tube and a Hopper 
shield (grey), a calibration unit attached to the tube. 
The Fully Coded (= arctan[(D M - D B )/L]), and the 
Partially Coded, (= arctan [ (Dm +Dd )/£]), Field of Views 
are indicated along with the telescope angular resolution 
(FWHM) 6 = arctan(#/L), where D D and Dm are the 
detector and mask dimensions, H is the width of the mask 
elements and L the distance mask-detector (« 3200 mm 
for ISGRI). 



by the convolution of the sky image S with M plus an un- 
modulated background array term B: D = S * M + B. If 
M has a correlation inverse array G such that M *G = 5- 
function, then we can reconstruct the sky by performing 
the following simple operation 

S' = D*G = S*M*G + B*G = S*5 + B*G = S + B*G 

and S' differs from S only by the B *G term. In the case 
the total mask M is derived from a cyclic replication of 
the same basic pattern and the background is given by 
a flat array £>, the term B * G is a constant and can be 
removed. Mask patterns with such properties, including 
the uniformly redundant arrays (URA), were found in the 
70s (Fenimore & Cannon 1978) and then succesfully em- 
ployed in X/7-ray telescopes onboard several high energy 
missions (Caroli et al. 1987, Paul et al. 1991). 

The IBIS coded mask imaging system (Fig. [Q in- 
cludes a replicated Modified URA (MURA)(Gottesman 
& Fenimore 1989) mask of tungsten elements (Fig. EJ) 
and 2 pixellated gamma-ray detector planes, both approx- 
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Fig. 2. The IBIS mask pattern of 95x95 elements (left) is 
formed by a replicated 53x53 MURA basic pattern, whose 
cyclic autocorrelation (right) is a 5 function. 



imately of the same size of the mask basic pattern: ISGRI, 
the low energy band (15 keV - 1 MeV) camera (Lebrun 
et al. 2003) and PICsIT, sensitive to photons between 
175 keV and 10 MeV, disposed about 10 cm below ISGRI 
(Di Cocco et al. 2003). The physical characteristics of the 
IBIS telescope define a FCFOV of 8° x 8° and a total 
FOV (FCFOV+PCFOV) of 19° x 19° at half sensitivity 
and of 29° x 29° at zero sensitivity. The nominal angular 
resolution (FWHM) is of 12'. ISGRI images are sampled in 
5' pixels while PICsIT images in 10' pixels. The MURAs 
are nearly-optimum masks and a correlation inverse is ob- 
tained by setting G = 2M - 1 (i.e. G = +1 for M = 1, 
G = — 1 for M = 0) apart from the central element which 
is set to 0. Simple correlation between such array and the 
detector plane array provides a sky image of the FCFOV 
where point sources appear as spikes of approximately the 
size (FWHM) of 1 projected mask element (12') with flat 
sidelobes (Fig. 2). The detailed shape of the System Point 
Spread Sunction (SPSF), the final response of the imag- 
ing system including the decoding process to the a point 
source, actually depends also on the instrument features 
(pixel size, detector deadzones, mask thickness, etc.) and 
on the decoding process (Gros et al. 2003). For this simple 
reconstruction (sum of transparent elements and subtrac- 
tion of opaque ones) and assuming Poissonian noise, the 
variance in each reconstructed sky pixel of the FCFOV is 
constant and simply given by V = G 2 *D = ^2 kl Dm, i- e - 
the total counts recorded by the detector. Therefore the 
source signal to noise is simply 



S/N = 



Cs 



y/Cs + Cb 



Reconstructed Source Counts 
\J Total Counts 



Sources outside the FCFOV but within the PCFOV still 
project part of the mask on the detector and their con- 
tribution can (and must) be reconstructed by properly 
extending the analysis in this region of the sky. In the 
PCFOV the SPSF does have secondary lobes (Fig. 01, the 
sensitivity decreases and the relative variance increases to- 
wards the edge of the field. In this part of the field also 
FCFOV sources will produce side lobes (see Sect. 2). 
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2. Sky image deconvolution 

The discrete deconvolution in FCFOV can be extended to 
the total (FC+PC) FOV by performing the correlation of 
the detector array D in a non cyclic form with the G array 
extended and padded with elements outside the mask. 
Since the number of correlated (transparent and opaque) 
elements in the PCFOV is not constant as in the FCFOV, 
the sums and subtractions for each sky position must be 
balanced and renormalized. This can be written by 

s = T,ki Gf +kij+l W k iDki _ Efcz Gr +k j+l W k iD k i 

>J J2kl G t+k,j+l W kl EfcZ G 7+k,j+l W kl 

where the decoding arrays are obtained from the G array 
by setting G+ = G for G > 0, G+ = for G < and 
G~ = -G for G < 0, G~ = for G > 0, and then are 
enlarged and padded with 0's outside the mask region. 
The sum is performed over all detector elements /c, / and 
z, j run over all sky pixels (both of FC and PCFOV). In 
the FCFOV we obtain the same result of the standard 
FCFOV correlation. To consider effects such as satellite 
drift corrections (see Goldwurm 1995), dead areas, noisy 
pixels or other instrumental effects which require different 
weighing of pixels values, a specific array W is used. For 
example W is set to for noisy pixels or detector dead 
areas and corresponding detector pixels are not included 
in the computation. The variance, which is not constant 
outside the FCFOV, is computed accordingly by 

_ E kl (Gt +k , j+l W k i) 2 D k i Y. kl {G- +Kj+l W kl ) 2 Dki 

13 (E kl Gt +kJ+l w kl y (E kl G~ +kJ+l w kl y 

since the cross-terms G + and G~ vanish. In the FCFOV 
the variance is approximately constant and equal to the 
total number of counts on the detector. The values com- 
puted by these formula are then renormalized to obtain 
the images of reconstructed counts per second from the 
whole detector for an on- axis source. Since mask elements 
are larger than detector pixels, the decoding arrays are 
derived as above after G has been sampled in detector 
pixels by projection and redistribution of its values on a 
pixel grid. For non integer sampling of elements in pixels 
(the IBIS case) G will assume continuous values from -1 to 
1. Different options can be adopted for the fine sampling of 
the decoding array (Fenimore & Cannon 1981, Goldwurm 
1995), the one we follow optimizes the point source signal 
to noise (while slightly spreading the SPSF peak). The 
SPSF, for an on-axis source and the IBIS/ISGRI con- 
figuration, obtained with the described deconvolution is 
shown in Fig. [31 Note the central peak and the flat level 
in the FCFOV, and the secondary lobes with the 8 main 
ghosts of the source peak at distances multiple of the mask 
basic pattern in the PCFOV. This procedure can be car- 
ried out with a fast algorithm by reducing previous formu- 
lae to a set of correlations computed by means of FFTs. 




Fig. 3. The System Point Spread Function over the com- 
plete FOV for the IBIS/ISGRI telescope. 



3. Non-uniformity and background correction 

In most cases and in particular at high energies the back- 
ground is not spatially uniform. Further modulation is 
introduced on the background and source terms by the 
intrinsic detector non-uniformity. The modulation is mag- 
nified by the decoding process (the term G * B in Sect. 
1) and strong systematic noise is generated in the de- 
convolved images (Laudet & Roques 1988). The result- 
ing systematic image structures have spatial frequencies 
similar to those of the modulation on the detector plane. 
For the IBIS detector, intrinsic non-uniformity is gener- 
ated by variations in efficiency of single detectors or asso- 
ciated electronics. Anticoincidence effects and scattering 
may also induce features at small spatial scales, of the or- 
der of the mask element size (Bird et al. 2003). The mod- 
ulation induced by the external background is normally a 
large scale structure dominated by the differences in effi- 
ciency of the veto system and of multiple event tagging. 
The background intensity may be variable on timescales 
from hours/days (activation in radiation belts, solar flares) 
to weeks/months (orbit circularization, solar cycle modu- 
lation) (Terrier et al. 2003). 

For the PICsIT layer, two more effects should be taken 
into account: the modulation at the border of the detec- 
tor semimodules (Di Cocco et al. 2003) caused by the 
onboard multiple reconstruction system, and the cosmic- 
rays induced events. The former is due to the fact that 
multiple events are reconstructed by analysing events per 
semimodule. When a multiple event is detected by two 
semimodules, it is considered as two single events. Since 
this occurs at the borders of the semimodules, there is an 
excess of single events in these zones, and a corresponding 
deficit of multiple counts in the same pixels. Obviously 
this effect becomes important as the energy increases. In 
Fig. 5 it is shown an example of a single event shadow- 
gram not yet corrected, where the border effects are clearly 
visible, and after the correction for background and non- 
uniformities (Natalucci et al. 2003). The second effect is 
due to spurious events caused by cosmic-rays interaction 
with Csl crystals (Segreto et al. 2003). This is relevant 
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Fig. 4. Results of reconstruction procedure applied to 
IBIS/ISGRI data of the Cygnus region: detector image af- 
ter background correction (up, left), associated efficiency 
image (up, right), decoded and cleaned sky image in the 
whole FOV (FC+PC) (down, left), associated variance 
(down, right). The algorithm detected, and identified 2 
sources: Cyg X-l at the center of the FOV and Cyg X-3 
in the PCFOV. 



only at low energies (< 250 keV) and results in a loss of 
sensitivity of about a factor of 4. It is possible to remove 
these fake events only if PICsIT is in photon-by-photon 
mode, while on the histograms it is possible to apply only 
an a-posteriori correction, by rescaling those pixels that, 
after the background subtraction, appear to be still noisy. 

If the shape of the background does not vary rapidly, 
regular observations of source empty fields can provide 
measures of the background spatial distribution which can 
be used to fiat-field the detector images prior to the decod- 
ing (contribution of weak sources is smoothed out by sum- 
ming images corresponding to different pointing aspects) 
(Bouchet et al. 2001). If U is the detector non-uniformity 
and B the background structure, then the recorded detec- 
tor image during an observation is D = (S*M)U + B and 
a basic correction can be performed by (D — bBEF)/Uc, 
where Bef is obtained from the empty field observations, 
Uq is an estimate of the detector non-uniformity from 
ground calibrations or Monte Carlo modelling and b is 
a normalization factor. If Uq is a good estimate of U and 
£?ef is close to B we obtain the needed correction. The 
normalization b is estimated either using the relative ex- 
posure times or the total number of counts in the images. 



4. The IBIS scientific data analysis 

The specific procedure for the IBIS data analysis starts 
from the data files obtained with the Integral Science Data 
Center (ISDC) (Courvoisier et al. 2003) preprocessing and 
performs a number of analysis steps which are hereafter 
described for the IBIS standard and photon-photon modes 
and the ISGRI and PICsIT data only. 

The ISDC preprocessing decodes the telemetry pack- 
ets, prepares the scientific and housekeping (HK) data 
with the proper satellite onboard times. This analysis task 
also computes the good time intervals by checking the 
satellite attitude, the telemetry gaps and the instrument 
modes and by monitoring the technological parameters. 
For the IBIS standard mode (Ubertini et al. 2003, Di 
Cocco et al. 2003) the following data sets are provided, 
for each elementary observation interval (science window) 
which corresponds to a constant pointing or slew satellite 
mode: 

- ISGRI event list: position, arrival time, pulse- height 
channel, rise-time. 

- PICsIT spectral-image histograms: 2 sets of 64 x 64 pixel 
images in 256 energy channels, one for the integrated sin- 
gle events and the other for the integrated multiple events. 

- PICsIT spectral timing histograms: count rates from the 
whole PICsIT camera integrated in a short time interval 
(0.97-500 ms, default: 2 ms) for up to 8 energy bands (de- 
fault: 4 bands). 

- Compton events list: ISGRI position, PICsIT position, 
ISGRI deposited energy, PICsIT deposited energy, ISGRI 
risetime, arrival time of events in coincidence between 
ISGRI and PICsIT. 

For the IBIS photon-photon mode the same ISGRI and 
Compton data are provided while PICsIT histograms are 
replaced by 

- PICsIT event list for single and double events: position, 
deposited energy, arrival time. 

The first tasks performed by the IBIS specific software are 
the computation of the livetimes of the single pixels, of 
the deadtimes and of the deposited energies for all events. 
Detector images are then built and corrected. Sky images 
are reconstructed for each pointing and then combined 
to obtain mosaic of sky images of the whole observation. 
Source spectra are extracted by performing image binning 
and correction on small energy bands and then compar- 
ing resulting images to source models. All data products 
are written in files with the standard FITS and OGIP for- 
mat in order to allow the use of standard high-level data 
analysis packages like FTOOLS, XSPEC and DS9. 

4.1. Pixels livetimes, deadtimes and energy corrections 

From the data of each single science window, prepared 
by the standard ISDC preprocessing, selected IBIS house- 
keeping (HK) data are first analyzed to prepare correction 
parameters for the scientific analysis. The HK tables re- 
porting the initial status (on/off) of detector pixels, their 
low energy thresholds and their status evolution during 
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the observations are decoded and the information orga- 
nized to allow computation of livetime of each single pixel 
during the science window. The deadtime and its evolu- 
tion is computed from the countrates reported in the HKs 
for the different IBIS datatypes, also including the time 
loss of random coincidences during veto rejection or cali- 
bration source tagging. 

The ISGRI event pulse height channel and risetime 
information in the science data list are then used to re- 
construct event deposited energy (in keV) by correct- 
ing the pulse height amplitude for the charge loss effect 
and for gains and offsets (Lebrun et al. 2003). The cor- 
rection is performed using look-up tables (LUT) derived 
from ground and inflight calibrations (Terrier et al. 2003). 
PICsIT energies are reconstructed on board using gain and 
offset LUT before histograms are accumulated. A further 
correction on the PICsIT single events can be performed 
on ground by using pixel gain and offsets values refined ac- 
cording to temperature variations (Malaguti et al. 2003). 

4.2. Image binning and uniformity-background 
correction 

For each pointing science window the event list or the his- 
tograms are then binned in detector images for the spec- 
ified energy bands and for a given risetime interval. The 
procedure computes an efficiency (or exposure) map com- 
bining informations from good time intervals, telemetry 
or data gaps, deadtimes, livetime of each pixel, and low 
energy thresholds. The map gives an exposure of each pix- 
els considering all these effects and relative to an effective 
exposure time of the pointing. 

ISGRI images are then enlarged to 130 x 134 pix- 
els and PICsIT images to 65 x 67 pixels to include the 
deadzones between the modules. Residual noisy pixels, i.e. 
pixels with large values (> 5 — 10 a) above the average 
are then identified and efficiencies of noisy and deadzone 
pixels are set to 0. The detector image pixels are divided 
by the non-zero-values of the efficiencies to renormalize 
their values. For efficiencies equal to the weighting array 
W is set to and the pixels are not included in the im- 
age or spectral reconstructions. The correction for detec- 
tor and background non-uniformity is then applied using 
the background images derived from empty field obser- 
vations and the detector uniformity maps (Terrier et al. 
2003, Natalucci et al 2003). For the spectral extraction 
however the correction is not perfomed at this level, back- 
ground and non-uniformity are accounted for directly in 
the fitting procedure. 

4.3. Iterative sky image reconstruction and cleaning 

From the corrected detector images the analysis procedure 
decodes the sky images using the algorithm described in 
Sect. 2 and iteratively searches for significant peaks in the 
image. The signal to noise levels (S/N) of detection can be 
set by the users, and the search can be preferencially per- 
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Fig. 5. Results of the background correction on a PICsIT 
single event shadowgram (252-329 keV) of the Crab re- 
gion. The raw image (left) shows presence of dead pixels 
and count excess on the module borders, not visible in the 
background corrected image (right). 

formed or even totally restricted to the sources included 
in an input catalogue. The search procedure is a delicate 
process due to the presence of the 8 ghosts (per source) 
and is sensitive to the choices of the S/N levels of search 
in particular if the background is not fully corrected. The 
first significant peak is detected and fitted with an ana- 
lytical approximation of the SPSF to finely determine the 
position of the source and its flux (Gros et al. 2003). The 
model of the shadowgram projected by the source at the 
derived position is computed, decoded, normalized to the 
observed excess and subtracted from the image. Such a 
source response model must take into account the mask el- 
ement thickness, the effects at the mask border, the differ- 
ent opacity of the mask support and of the passive shield, 
and the transparent zones in or between these elements. 
However the model includes only the geometrical absorp- 
tion effects. Proper account of scattering and energy re- 
distribution is included in the energy response. The same 
routine which computes the model of source shadowgram 
on the detector is used in the spectral extraction. The 
source parameters (position, flux, and signal to noise) are 
stored in an output file and the procedure uses the source 
catalogue to identify the source. A second excess is then 
searched and cleaned and the procedure continues till all 
sources or excesses are identified and cleaned. The image 
of the main peak of each source is finally restored in the 
sky image while the secondary lobes are excluded. Units 
of intensity images are in reconstructed counts per second 
for an on-axis source, and a variance image is computed 
in the proper units. 

4.4. Image mosaic and final image analysis 

The reconstructed and cleaned sky images of each pointing 
science window are then rotated, projected over a refer- 
ence image and summed after a properly weighting with 
their variance and exposure time. There are two ways to 
implement the rotation of an image. Either the pixel value 
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Fig. 6. Sky image mosaic of the Cygnus region in the 15- 
40 keV band composed with IBIS/ISGRI data from 130 
pointings for a total exposure of 145 ks. Three sources are 
detected at high significance level: Cyg X-l, Cyg X-3 and 
EXO 2030+375. 



is redistributed around the reference pixels on which is 
projected or its value is entirely attributed to the pixel 
whose center is closer to the projected center of the in- 
put pixel. In the first case the source is slightly smoothed 
and enlarged but the source centroid is well reconstructed, 
while in the second case discrete effects can produce a bias 
in the source centroid but the intensity will be better re- 
constructed and the source less spread. 

After the mosaic step is completed the routine restart 
the search and analysis of point sources in the mosaicked 
sky image, again using a catalogue for the identification. 
Results of fluxes and positions are reported in output. 
Since the point source location error is inversely propor- 
tional to the source signal to noise (Gros et al. 2003), loca- 
tion of weak sources is better carried out on mosaicked im- 
ages. Note that the search for significant excesses (both for 
single pointing and for the mosaic) must be performed tak- 
ing into account that these are correlation images and con- 
sidering the number of independent trials that are made by 
testing all pixels of the reconstructed sky image. The criti- 
cal level at which an unknown excess is significant must be 
increased from the standard 3 a value to typically 5 - 6 a 
(see e.g. Caroli et al. 1987). In the case of residual (back- 
ground) systematic modulation, the detection level must 
be further increased (for example by a factor given by the 
ratio of measured to computed image standard deviation). 



4.5. Source spectra and light curve extractions 

Once the positions of the active sources of the field are 
known, their fluxes and count spectra are extracted for 
each pointing science window in predefined energy bins. 



Detector images and their efficiencies are built for these 
energy bins, enlarged and then searched for residual noisy 
pixels. For each source and each energy band a shadow- 
gram model is then computed using the detailed modelling 
algorithm described above. If P kl (AEh) is the model for 
the source n and the energy bin AEh, the total sky plus 
background model T at the pixel fc, I is defined as 



T%(AE h ) = J2 tiP^E h )E h kl U% + b h B kl (AE h )E- 



J ki 



where the sum is extended to all "active" sources of the 
region and where the background image B is derived from 
the empty field data, the Uq is the estimated detector non- 
uniformity and E is the efficiency image. This model is fit- 
ted to the enlarged and noisy-pixel-corrected detector im- 
ages Dk,i by using the maximum likelihood technique with 
Poissonian distribution. The optimization gives the mul- 
tiplicative factors f% and b h for each energy band which 
(after proper normalizations and division by the exposure 
time), provide the source and background count rate spec- 
tra. It is important to perform the simultanous fit for the 
models of all the active sources of the field to avoid con- 
tamination by other sources of the interested spectrum. 
Spectra of the same source collected in different point- 
ing science windows can then be summed to obtain an 
average spectrum during an observation. The spectra are 
compared to physical spectral models convolved with the 
energy response function of the instrument (Laurent et 
al. 2003) in order to evaluate source physical parameters. 
Finally, a similar procedure, where the images are binned 
in large energy bands and small time intervals, provides 
source light curves on time bins shorter than science win- 
dows for the predefined active sources in the field of view. 
For PICsIT operating in standard mode, light curves of 
the count rates registered by the whole detector can be 
generated directly from the spectral timing data and used 
to search for characteristic variability signatures like pul- 
sations or bursts. 

5. Final remarks 

The software has been tested using inflight data in several 
conditions and the performance is satisfactory consider- 
ing the early phase of the mission. Background corrections 
have been already implemented for the analysis of PICsIT 
data while not yet fully used for the ISGRI analysis (but 
see Terrier et al. 2003). Fig. shows the correction oper- 
ated on the Crab nebula PICsIT images as described in 
Sect. 3. In Fig. the resulting sky image mosaic obtained 
from the IBIS scientific pipeline applied to about 130 el- 
ementary pointings on Cygnus X-l, illustrate the qual- 
ity of the scientific analysis presently available. Images 
are well reconstructed and cleaned and sources well de- 
tected and finely positioned, with typical absolute error 
radii (90% confidence level) of less V for strong source 
(S/N > 30-40) and less than 3' for moderate or weak ones 
(S/N « 10-40) (Gros et al. 2003). Fig. Q shows a recon- 
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Fig. 7. IBIS/PICsIT reconstructed image of the Crab in 
the energy band 252-329 keV obtained from the IBIS sci- 
entific pipeline applied to 55 science windows for a total 
of 132 ks of exposure. 
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Fig. 8. IBIS/ISGRI reconstructed count spectrum of the 
Crab from data of a 90 ks exposure pointed observation. 
The best fit model of a power law and ratios are also 
shown. The ratios show presence of 10% systematic erros 
in particular in the 20-100 keV range. 



structed PICsIT image of the Crab, a very bright con- 
stant and pointlike source for IBIS, often used as calibra- 
tion source for high energy instruments. Fig. [HI reports a 
reconstructed Crab count spectrum obtained with the de- 
scribed scientific analysis procedures from the ISGRI data 
of an observation pointed on the source and performed in 
staring mode (constant attitude asset during the observa- 
tion, see Winkler et al. 2003). The spectrum is compared 
to a model of a simple power law using the instrument 
energy response files. The best fit in the range 20-700 keV 
is found for a slope of 2.05 which is compatible with the 
high energy spectrum of the Crab. To have acceptable x 2 > 
systematic errors of about 10% have to be included. 

These examples give an idea of the performance of 
the analysis s/w and of the calibration and response files 
presently available. The procedures described here have 
been integrated in the ISDC analysis system and consti- 
tute the core of the IBIS data analysis software of the 
Quick Look and Standard Analysis run at the ISDC and 
of the INTEGRAL Data Analysis System package. The 
IBIS analysis software and associate calibration files are 
constantly improved and integrated in the following ver- 
sions of the ISDC system. The present set has already 
been used to derive a number of interesting results from 
the IBIS data of the high energy sources observed with 
INTEGRAL (see this volume). 
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